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XO-2b: a Prograde Planet with a Negligible Eccentricity, 
and an Additional Radial Velocity Variation* 
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Abstract 

We present precise radial velocities of XO-2 taken with the Subaru HDS, covering two transits of X0-2b 
with an interval of nearly two years. The data suggest that the orbital eccentricity of X0-2b is consistent 
with zero within 2a (e = 0.045 ±0.024) and the orbit of X0-2b is prograde (the sky-projected spin-orbit 
alignment angle A — 10° ± 72°). The poor constraint of A is due to a small impact parameter (the orbital 
inclination of X0-2b is almost 90°). The data also provide an improved estimate of the mass of X0-2b as 
0.62 ±0.02 Afjup. We also find a long-term radial velocity variation in this system. Further radial velocity 
measurements are necessary to specify the cause of this additional variation. 
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1. Introduction 

Planetary orbits in binary systems provide useful clues 
to learn planetary migration mechanisms in binary sys- 
tems. Planets in binary systems have a chance to take a 
different path of orbital migration than planets in single 
star systems. That is, due to the presence of a binary com- 
panion, an inner planet can evolve into a highly eccentric 
and highly tilted, or even retrograde orbit, by the mech- 
anism known as the Kozai migration (Kozai 1962; Wu 
& Murray 2003; Fabrycky & Trcmaine 2007; Lithwick 
& Naoz 2011; Katz et al. 2011). Observable evidences 
of the Kozai migration are large eccentricity and large 
spin-orbit misalignment of an inner planet. One can mea- 
sure the orbital eccentricity via radial velocity (RV) mea- 
surements, and one can also learn the spin-orbit align- 
ment angle at least in sky-projection via the Rossiter- 
McLaughlin effect (hereafter the RM effect: Rossiter 1924, 
McLaughlin 1924; Ohta et al. 2005; Hirano et al. 2010; 
Hirano et al. 2011a) in transiting systems. 

One interesting example of a transiting planet in a wide 
binary system is HD 80606b, which has a large eccentricity 
(e = 0.93, Naef et al. 2001) and a large sky-projected spin- 
orbit alignment angle (A = 42°, Moutou et al. 2009; Pont 
et al. 2009; Winn et al. 2009b; Hebrard et al. 2010). The 
orbit of HD 80606b has been well explained by a scenario 
of the Kozai migration due to the presence of the binary 
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companion HD 80607 (Wu & Murray 2003). Another 
interesting example is WASP-8b, which has an eccen- 
tric and retrograde orbit with a wide binary companion 
(Queloz et al. 2010), although a Kozai migration scenario 
for WASP-8b has not been well examined. Other tran- 
siting planets in binary systems, however, show neither 
large eccentricity nor large spin-orbit misalignment (e.g., 
HD189733b: Bakos et al. 2006; Winn et al. 2006, TrES- 
4b: Mandushev et al. 2007; Narita et al. 2010b). This 
fact can be understood by considering that the Kozai mi- 
gration needs to meet some stringent conditions (see e.g., 
Innanen et al. 1997; Wu et al. 2007; Narita et al. 2010c). 
To learn the occurrence frequency of the Kozai migration 
for planets in binary systems, it is important to increase 
the number of observations for such transiting systems. 

We here report on the orbit of the transiting planet 
X0-2b, which was discovered by Burke et al. (2007) in 
the course of the XO survey. The host star XO-2 is a 
relatively faint {V = 11.18) KOV star (the stellar mass 
0.98±0.02Ms and the stellar radius 0.97±0.02i?s) at -150 
pc from the Sun (Burke et al. 2007). XO-2 has the twin 
(the same spectral type) stellar companion X0-2S at 31" 
(— 4600AU) separation. We note that the Kozai migration 
for X0-2b could occur if the eccentricity of the binary 
orbit (currently unknown) is very large (e > 0.9), when 
the timescale of the Kozai migration is shorter than the 
timescale of a perturbation due to General Relativity for 
X0-2b (Wu et al. 2007). Burke et al. (2007) reported the 
mass, radius, orbital period of X0-2b as 0.57±0.06 Mjup, 
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tively. The eccentricity of X0-2b was assumed to be zero 
in the discovery paper, and no measurement of the RM 
effect in this system has been reported. In this letter, we 
present RVs of XO-2 spanning about two years and cov- 
ering two full transits of X0-2b allowing us to model the 
RM effect. Our RV data suggest a small eccentricity and 
a likely spin-orbit alignment of X0-2b, concluding no sup- 
porting evidence of the Kozai migration. In addition, we 
find a long-term RV trend for XO-2. This RV acceleration 
may suggest a third body in this system. 

The rest of this letter is organized as follows. We 
present details of our observations and reductions of RV 
data in section 2. We describe our model to fit the ob- 
served data in section 3. We present our results on the 
orbit of XO-2 and discussions on possible migration mech- 
anisms of X0-2b in section 4. Finally, we summarize the 
findings of this letter in section 5. 

2. Observations and Data Reductions 

We observed XO-2 with the standard I2a setup of the 
High Dispersion Spectrograph (HDS: Noguchi et al. 2002) 
onboard the Subaru 8.2m telescope on Mauna Kea. An 
iodine gas absorption cell was inserted for RV measure- 
ments. We measured RVs around two full transits of XO- 
2b on UT 2008 March 9 and UT 2009 November 24. We 
also gathered out-of-transit RVs spanning about two years 
from 2008 to 2010. Slit width was either 0."4 {R = 90000) 
or 0'.'6 {R = 60000) depending on observing conditions. 
Exposure times were 600 s for around-transit phase and 
900 s for out-of-transit phase. Finally, we obtained a high 
signal-to-noise ratio (SNR over 200) and high spectral 
resolution (R = 90000) template spectrum on UT 2011 
February 15 with an exposure time of 1800 s. 

We process the observed frames with standard IRAF 
procedures and extract ID spectra. Relative RVs and un- 
certainties are computed by the algorithm of Butler et al. 
(1996) and Sato et al. (2002). We estimate the uncer- 
tainty of each RV based on the scatter of RV solutions for 
~4 A segments of each spectrum as described in Narita 
et al. (2007). The RVs and internal errors are summa- 
rized in table 1, and plotted in the top panel of figure 1. 
Note that we use HJDutc (Heliocentric Julian Date in 
the Coordinated Universal Time) as a time standard. 

Since we do not have good photometric transit data, we 
additionally incorporate five published photometric light 
curves taken with Keplercam on the 1.2m telescope at the 
Fred Lawrence Whipple Observatory (FLWO) on Mount 
Hopkins, Arizona (Fernandez et al. 2009). Those tran- 
sits were observed in the Sloan z' band. Following the 
procedure by Fernandez et al. (2009) and Narita et al. 
(2010a), we rescale the photometric uncertainties to ac- 
count for time-correlated noise (so-called red noise: see 
e.g., Pont et al. 2006) by calculating a red noise factor 
/? = crN.obs/cN. ideal for various N (corresponding to 10-20 
min), and multiply the photometric uncertainties of each 
dataset by the maximum value of /3. We note that we find 
slightly higher values of /3 than Fernandez et al. (2009), 
and thus our photometric uncertainties are conservative. 



3. Model 

First we model the RVs assuming a single orbiting 
planet. We then find RVs around the second transit are 
vertically off by about 15 m s~^ from those around the 
first transit (see the middle panel of figure 2). We also find 
that residuals from the single planet model has a linear 
trend in time (see the middle panel of figure 1). We thus 
model our data assuming a single orbiting planet with 
a linear RV acceleration. The five FLWO transit light 
curves arc simultaneously fitted assuming the quadratic 
limb-darkening law. We note that we fix one coefficient 
(iti) of the limb-darkening parameters to ui = 0.25, which 
is the central value derived by Fernandez et al. (2009). 
This treatment is useful so as to avoid underestimate of 
uncertainties for fitted parameters (see Southworth 2008). 

Consequently, our model has 16 free parameters: the 
RV semiamplitude K, the eccentricity e, the argument of 
periastron zu, the sky-projected stellar rotational velocity 
T^sin/s, the sky-projected spin-orbit alignment angle A. 
the offset RV for the Subaru data 7, the RV acceleration 
7, the planet-star radii ratio Rp/Rg, the orbital inclina- 
tion i, the semi- major axis in units of the stellar radius 
a/Rs, one of the limb-darkening coefficients U2, and the 
five mid-transit times of the FLWO data Tc{E). We fix 
P = 2.6158640 days and Tc{0) = 2454466.88467 in HJDutc 
which are the values reported by Fernandez et al. (2009). 
Uncertainties of mid-transit times for the two spectro- 
scopic transit observations (UT 2008 March 9 and UT 
2009 November 24) are well within the shortest exposure 
time of RVs (600 s) and have little effect on results. 

The statistic for the joint fit is 
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where /i,obs, Vj.ohs, Cj, and aj are the observed relative 
fluxes, RVs, and their uncertainties. We note that we 
do not add RV jitter in aj, since reduced for the RV 
data is below unity as shown later (see table 2). The 
modeled fluxes (/i,caic) are computed by the formula given 
by Ohta et al. (2009), and the modeled RVs (wj,modei) are 
calculated as Wcaic = WKopicr + wrm +7^ + 7; where WKopier 
is the Keplerian motion, wrm is the RM effect, and jt 
represents the long-term RV trend. We compute the RM 
effect URM by the accurate analytic formula for a KOV 
star described in Hirano et al. (2011a). In the formula, we 
use ui = 0.714 and U2 = 0.114 for the band of the iodine 
absorption lines, based on the table of Claret (2004). We 
then determine optimal parameter values by minimizing 
the statistic using the AMOEBA algorithm (Press et al. 
1992). Uncertainties of free parameters are estimated by 
the criterion Ax^ = 1.0. 

4. Results and Discussions 

The bcst-flt parameters and errors are summarized in 
table 2. Photometric parameters related with only the 
FLWO data {Rp/Rs, 112, and Tc) are essentially consis- 
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Fig. 1. Top panel: RVs of XO-2 observed with the Subaru 
HDS. Middel panel: Residuals of RVs from the best-fit model 
without subtracting the long-term RV trend. Bottom panel: 
Same as the middle panel but with subtracting the RV trend. 
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Fig. 2. Top panel: Phased RVs around transits after sub- 
tracting the RV trend. The circle and square symbols repre- 
sent RVs on UT 2008 March 9 and UT 2009 November 24, 
respectively. Middel panel: Residuals of RVs from the best-fit 
model without subtracting the long-term RV trend. Bottom 
panel: Same as the middle panel but with subtracting the RV 
trend. 



> 




-0.4 -0.2 



Orbital phase 



Fig. 3. The same as figure 2, but for a whole orbital phase. 



tent with the results of Fernandez et al. (2009) and are 
not shown. Figure 1 plots the Subaru RVs (the top 
panel) , residuals from the best- fit model without (the mid- 
dle panel) and with (the bottom panel) a subtraction of 
the RV trend jt. The for 59 RVs is 48.8 and the rms 
of RVs is 3.08 m s~^. This validates our treatment of 
stellar jitter in the RV errors. We find that the RV semi- 
amplitude K (92.2 ± 1.7 m s~^) is about 8% larger than 
than the value (85 ± 8 m s^^) reported by Burke et al. 
(2007) with a smaller uncertainty. This provides an im- 
proved estimate of the mass of X0-2b as 0.62 ±0.02 Afjup. 
Wc also find the eccentricity of X0-2b is small and consis- 
tent with zero within about 2a. We note that the results 
and errors are essentially unchanged even if we jointly fit 
the RV data by Burke et al. (2007). 

Figure 2 shows phased RVs around the transit phase 
with the best-fit model. The RV trend is subtracted in the 
top and bottom panels, but not subtracted in the middle 
panel. Figure 3 is the same as figure 2, but for a whole 
orbital phase. The shape of the RM effect clearly suggests 
a prograde orbit of XO-2. However, the derived constraint 
on A is very poor, namely A = 10° ± 72°. It is because the 
orbital inclination i of X0-2b is near 90°. Figure 4 plots 
a Ax^ contour map in (A, ysinJg) space. Our best-fit 
value of Fsin/g based on the RM effect is 1.45 km s~^, 
which is in good agreement with a value 1.4 ±0.3 km s~^ 
derived by Burke et al. (2007) based on a spectroscopic 
line analysis. Schlaufman (2010) independently estimated 
Vsin/s of XO-2 as 1.69 ±0.2 km based on its mass and 
age, and concluded that the XO-2 system is likely to be 
aligned. The conclusion is apparently consistent with our 
result. Thus although a larger value of Fsin/g is allowed 
by the RM effect, such a larger Fsin/g is unlikely when 
considering the other analyses. As a test case, if we add 
a prior constraint 1.4 ±0.3 km s~^ to the statistic, we 
find A = 9°t^^I . Thus the orbit of X0-2b is unlikely to be 
highly tilted, but an exact value of A is unfortunately in- 
determinate. As a result, we conclude a small eccentricity 
and a likely spin-orbit alignment for the orbit of X0-2b. 
As is known today, RM measurements have shown that 
about two thirds of close-in giant planets are spin-orbit 
aligned and prograde, whereas the last third is strongly 
misaligned or even retrograde (e.g., Hebrard et al. 2011). 
Our result suggests that X0-2b is apparently in the first 
category. 

We determine the RV acceleration of XO-2 as 7 = 
7.51 ± 0.58 m s^^ yr~^. The trend is small but signif- 
icant in the two years. If we do not include 7 as a 
free parameter, the value for the RVs changes from 
48.8 to 231.9. Assuming a long-term linear trend, this 
RV acceleration can be explained by a hypothcrical third 
body (index "c") whose mass and orbit follow a relation 
Mc sin ic/a^ =0.042 ±0.003 Mjup AU"^ The orbital pe- 
riod of the third body can be longer than ^8yr (^ 4AU) 
since the RV trend is apparently linear spanning about 
2 years. For example, a 1 Afjup planet at 5 AU like the 
Jupiter in our Solar System meets the relation. We note 
that the binary companion X0-2S does not meet this re- 
lation. There is another possibility that the timescale of 
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the additional RV variation is a few hundred days, due 
to the lack of our observations between the two transits. 
Additional RV measurements with the Subaru HDS will 
allow us to discriminate the timescale of the RV variation. 

Other possibilities of the trend include a systematic 
RV variation of the Subaru HDS, starspots or magnetic 
cycles of XO-2. The first scenario is unlikely, because 
RVs of the Subaru HDS for RV standard stars 
ble within a few m s~^ (Harakawa et al. 2010; Sato et al. 
2009). The starspot scenario is also unlikely, because no 
significant RV variation has been observed for XO-2 in 
a timescale corresponding to the rotation period of the 
host star (~35 day: Matsumura et al. 2010), although 
RV variations caused by spots should have been maxi- 
mum. At this point, we cannot exclude the possibility 
of magnetic cycles of XO-2 recently examined by Lovis 
et al. (2011), although any strong magnetic activity has 
not been reported for this system. To specify the cause of 
the RV variation and make a decisive conclusion, further 
RV monitoring would be very important in this system. 

There is still a possibility that the mutual inclination 
between the orbital axes of the transiting planet X0-2b 
and the binary companion X0-2S was once larger than the 
threshold value of the Kozai mechanism. In addition, if 
the eccentricity of the binary orbit is larger than ^-^0.9, the 
timescale of the Kozai migration might have been shorter 
than that of General Relativity for X0-2b at the birth- 
place (~2Gyr based on Wu & Murray 2003). However, the 
presence of the hypotherical third body completely dic- 
tates that X0-2b cannot be migrated through the Kozai 
migration caused by X0-2S, since the gravitational per- 
turbation timescale due to the third body (comparable to 
the orbital period of the third body) is much shorter than 
the Kozai migration timescale (see e.g., Innanen et al. 
1997; Narita et al. 2010c). 

On the other hand, Winn et al. (2010) pointed out that 
a hot Jupiter around a cool (less than '^6250K) star can 
lead to a re-alignment of its spin-orbit alignement angle 
due to the tidal force on a convective surface layer of the 
cool host star. XO-2 is in the category of a cool star 
(^5340K: Burke et al. 2007) and the spin axis of the host 
star might have been re-aligned. Thus one possible migra- 
tion scenario is that X0-2b has migrated through planet- 
planet scattring (e.g., Rasio & Ford 1996; Nagasawa et al. 
2008; Chatterjee et al. 2008), and subsequently the ec- 
centricity and the spin-orbit alignment angle have been 
damped. In that case, the hypothetical third body may be 
the counterpart of the planet-planet scattering. Another 
possible migration mechanism is the standard disk-planet 
interaction mechanism (e.g., Lin et al. 1996; Ida & Lin 
2004). Although we cannot specify the migration mech- 
anism of X0-2b at this point in time, further measure- 
ments of transit timings of X0-2b would be useful to con- 
strain another low-mass body in mean motion resonance, 
which can be an evidence of the disk-planet migration 
(Fernandez et al. 2009). Any detection or constraint of 
such a mean motion resonance body by transit timings, 
as well as further RV measurements, would allow further 
discussions on the migration history of this system. 
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Fig. 4. A contour map in X-Vsinis space. Tiie solid lines 
represent contour for Ax^ = 1-0. 

5. Summary 

We have monitored RVs of XO-2 spanning about two 
years with the Subaru HDS. We find a small eccentric- 
ity and a likely spin-orbit alignment for the orbit of the 
transiting planet X0-2b, and we also detect a long-term 
RV acceleration of the host star. Based on the observed 
properties of this system, we constrain the Kozai migra- 
tion scenario of X0-2b. This illustrates that a presence of 
a wide binary companion does not necessarily suggest the 
Kozai migration (see also Narita et al. 2009; Winn et al. 
2009a; Narita et al. 2010c, for the migration mechanism 
of a retrograde planet HAT-P-7b). To constrain the oc- 
currence frequency of the Kozai migration for planets in 
binary systems, combinations of observations of (1) the or- 
bital eccentricity by RV measurements, (2) the spin-orbit 
alignment angle by the RM effect, and (3) the binarity (or 
multiplicity) of the system by high-contrast direct imag- 
ing, are of special importance. Accumulating those obser- 
vations for transiting planetary systems would allow us to 
learn more about planetary migration mechanisms in the 
future. 
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Table 1. RV data taken with the Subaru HDS. 



SUbaru HDS RVs 
Time [HJD.UTC] Value [m s^^] Error [m s~^] 



2454534.71592 


0.00 


3.53 


2454534.72354 


5.83 


3.09 


2454534.73116 


-0.19 


3.03 


2454534.74713 


2.47 


3.28 


2454534.75474 


-4.44 


2.92 


2454534.76237 


-9.72 


2.68 


2454534.76999 


-8.87 


2.88 


2454534.77761 


-8.62 


2.73 


2454534.78523 


-6.95 


3.11 


2454534.79284 


-8.66 


2.85 


2454534.80046 


-12.98 


2.82 


2454534.80807 


-19.25 


3.35 


2454534.81570 


-14.86 


2.92 


2454534.82333 


-17.81 


3.20 


2454534.83095 


-17.81 


2.76 


2454534.83856 


-21.79 


3.16 


2454534.84618 


-18.34 


3.07 


2454534.85380 


-19.09 


3.19 


2454534.86141 


-15.27 


2.89 


2454534.86904 


-16.67 


3.85 


2454534.87665 


-24.82 


3.32 


2454534.88427 


-23.42 


3.00 


2454534.89188 


-26.77 


3.12 


2454534.89949 


-37.70 


3.60 


2454534.90711 


-39.15 


3.17 


2454534.91472 


-44.97 


3.35 


2454534.92234 


-50.05 


3.42 


2454534.92996 


-51.32 


3.43 


2454534.93759 


-50.14 


3.05 


2454534.94520 


-51.64 


2.84 


2454534.95282 


-43.37 


3.53 


2454534.96043 


-48.30 


3.25 


2454534.96805 


-48.68 


3.22 


2454616.73649 


-120.61 


3.04 


2454616.75121 


-124.63 


3.02 


2455160.01806 


-6.98 


2.87 


2455160.02580 


-10.17 


3.44 


2455160.03354 


-2.75 


3.50 


2455160.04127 


-0.66 


3.28 


2455160.04900 


-3.77 


3.53 


2455160.05675 


3.62 


3.52 


2455160.06449 


-7.39 


3.23 


2455160.07224 


-8.94 


2.90 


2455160.07999 


-16.07 


2.96 


2455160.08773 


-16.25 


3.23 


2455160.09570 


-28.54 


3.52 


2455160.10344 


-34.43 


2.95 


2455160.11118 


-34.92 


2.78 


2455160.11893 


-38.27 


3.51 


2455160.12668 


-35.58 


2.85 


2455160.13442 


-39.42 


2.75 


2455160.14217 


-32.60 


2.98 


2455160.14991 


-32.10 


2.83 


2455210.99190 


-45.61 


3.11 


2455211.00660 


-39.25 


2.88 


2455211.02128 


-44.23 


3.70 


2455211.75316 


63.93 


4.56 


2455211.76263 


71.23 


4.44 


2455232.09437 


3.96 


3.67 



Table 2. Best-fit values and errors of the parameters. 



Free Parameter 


Value Error 


K [m s-^] 
e 

vu [°] 

V smlg [km s J 
A[°] 

Ml 
a/Rs 

r — 1 1 

7 [m s J 

7 [m s^"'" yr^^] 


92.2 ±1.7 
0.045 ±0.024 

270 _7 

1 Ar- 4-2 73 

1.45 toll 
10 ±72 
88.7 tli 
8.43 

-34.3 ±1.1 
7.51 ±0.58 


Derived Parameter 


Value Error 


Mp [Mjup] 

McSinic/a? [Mj^p AU"^] 

for 59 RVs 
rms [m s~^] 


0.62 ±0.02 
0.042 ±0.003 
48.8 
3.08 



* All data are presented in the electric table. 



